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ABSTRACT 

We searched the literature to obtain a complete list of known Class sources. A list of 95 confirmed 
or candidate objects was compiled. To the best of our knowledge, all published broad-band observations 
from 1 fim to 3.5 mm have been collected and are assembled in a catalogue. These data were used to 
determine physical properties (Tboi, Lboi, Lgmm/Lbob Mcnv) and for a uniform classification. 50 sources 
possess sufficient observational data and are classified as Class or Class 0/1 objects. The source proper- 
ties are compared with different evolutionary models to infer ages and masses, and their correlations are 
investigated. About 25% of the sources are found to be in a quiet accretion phase or possess a significantly 
different time evolution of the accretion rate than the average. In Taurus, with its isolated star formation 
mode, this seems especially to be the case. 

Subject headings: Catalogs - Stars: evolution - Stars: formation - Infrared: stars 



1. Introduction 

Protostars begin their life soon after the collapse of 
their parental cloud core begins. Such stellar embryos 
have masses of less than 10^^ Mq (Larson (2003)) and 
are deeply embedded in a dense and massive envelope 
of gas and dust. The central object grows in mass by 
accreting material from this envelope either by direct 
infall or via a flattened circumstellar disc. As long as 
the central object is less massive than the surround- 
ing envelope it is called a Class object. Otherwise it 
is a Class I protostar. Since this mass ratio is not di- 
rectly observable, other criteria have to be adopted for 
a classification. Andre et al. (2000) defined three ob- 
servational properties needed to classify an object as 
a Class source: (i) an internal heating source (com- 
pact centimeter radio continuum) or bipolar outflow 
as an indication of a central young stellar object; (ii) 
extended and centrally peaked sub-millimeter contin- 
uum emission, indicating the spheroidal envelope; (iii) 
a high ratio of sub-millimeter to bolometric luminos- 
ity (Lsmm/Lboi > 0.005), with Lsmm measured long- 
ward of 350 /zm. It is shown in Andre et al. (1993) 
that this last criterion is equivalent to the mass ratio of 
the envelope and the central star being larger than one. 



Other authors use the bolometric temperature (Tboi) 
instead of Lsinm/Lboi (e.g. Chen et al. (1995a; 1997) 
use Tboi < 70 K). 

The Class stage can be characterised as the main 
mass accretion phase, in which the forming star gains 
the bulk of its final mass. This process is not yet well 
understood. How does the mass accretion rate depend 
on the age and the properties of the cloud, and is it 
governed by turbulence or ambipolar diffusion? What 
is the lifetime of the Class sources? How are the 
source properties connected to the molecular outflow, 
that is inevitably driven by these objects? The answers 
to these questions will finally help us to understand the 
feedback of the star formation process to the surround- 
ing molecular cloud and the initial mass function. 

To seriously address these questions we need a 
large, homogeneously classified sample of these very 
young protostellar objects. Several samples have been 
compiled in the recent years (e.g. Chen et al. (1995a), 
Andre et al. (2000), Shirley et al. (2000), Motte & 
Andre (2001)). All these works have specific prob- 
lems: (1) they concentrate only on a limited number 
of sources; (2) different criterias for classification of 
Class objects are used; (3) only parts of the available 
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observational data are considered. 

The first step to obtain a large, homogeneously clas- 
sified sample of Class sources is a complete search 
of the Uterature for such objects. Here we present an 
as far as possible complete list of the youngest known 
protostars to date. We will classify them by means 
of their spectral energy distribution (SED) using all 
(to the best of our knowledge) published broad-band 
photometric data. We derive basic properties of these 
sources (Tboh Lboi)- The sub-mm data are used to es- 
timate envelope masses (Mcnv) and sub-mm slopes of 
the SED (/3). By comparing Tbob Lboi, and Mcnv with 
evolutionary models we determine ages and masses of 
these sources and investigate limitations of the current 
observational data and evolutionary models. This is of 
particular interest for forthcoming powerful telescopes 
like Spitzer and ALMA, which will enable us to im- 
prove our knowledge of these sources significantly. 
Given this expected growth in observational data in 
the near future, the present catalogue might be quickly 
outdated. Hence, all the data will be available on a 
webpage' and updated regularly. 

In Sect. 2 we introduce the source sample. The 
method of our data analysis is then put forward in 
Sect. 3, including a description of the source classifi- 
cation in Sect. 3.2. Finally the results and implications 
are discussed in Section 4. 

2. Source sample and broad band photometry 

To obtain an as far as possible complete list of 
known Class objects we combined the samples of 
candidates and sources mentioned in various publica- 
tions. In particular we combined: (1) sources with 
Tboi < 100 K from Chen et al. (1995a); (2) objects 
with Tdust < 100 K from Hurt et al. (1996); (3) Class 
sources listed in Bontemps et al. (1996); (4) objects 
of type CLO in Saraceno et al. (1996); (5) candidates 
with Tboi < 100 K from Chen et al. (1997); the source 
B 35 could not be identified in SIMBAD; (6) Objects 
with Lboi/Lsmm <200 in Chini et al. (1997a) with 
VLA detected counterparts in Reipurth et al. (1999); 
(7) the Class sources in CB 68 and CB 232 men- 
tioned in Huard et al. (1999); (8) objects classified as 
Class in Park et al. (1999); (9) all objects listed in 
Andre et al. (2000); (10) Class sources in Shirley 
et al. (2000); (11) sources indicated as Class or 
Class 0/1 in Motte & Andre (2001); (12) objects with 

^http://www.dias.ie/protostars 



LFiR/Lsmm < 200 in Chini et al. (2001); (13) candi- 
dates with Tboi < lOOK in Lehtinen et al. (2001); (14) 
Class sources from Visser et al. (2002); (15) objects 
with Tboi < 100 K Hsted in Young et al. (2003); (16) 
Class objects and candidates investigated in Froe- 
brich et al. (2003); (17) candidates and confirmed 
Class sources from Rengel et al. (2004); (18) the 
possible Class sources NGC7129FIRS2 (Eiroa et 
al. (1998)), NGC 2068 LBS 17 (Gibb & Little (2000)), 
IC348MMS (Eisloffel et al. (2003)), NGC 7538 S 
(Sandell et al. (2003)), and MonOB IIRAS 1 2 S 1 
(Wolf-Chase et al. (2003)). All investigated sources 
with their positions, adapted distances, and references 
are listed in Appendix A in Table 3. This table will be 
only available in electronic form and lists the object 
name used in this paper and other common names of 
the sources. 

Sources in the combined sample were identified by 
their position in the SIMBAD database in order to 
avoid confusion due to unclear or conflicting nomen- 
clature. A search in the related publications given 
by SIMBAD for the sources was performed and to 
our best knowledge all available broad-band photo- 
metric data were extracted. We restricted the search 
for data ranging from the near infrared (NIR) to wave- 
lengths shorter than 3.5 mm. At longer wavelengths 
the emission of the envelope might be a combination 
of dust continuum and free-free emission from the cen- 
tral source and may influence the analysis of the enve- 
lope properties. NIR and mid-infrared detections with 
brightness given in magnitudes were converted into 
Jansky using flux zero-points of Wamsteker (1981). In 
case of non-existing NIR detections in the literature 
we searched the 2MASS catalogue. If the source was 
not detected by 2MASS, we determined an upper limit 
(the faintest flux in the 2MASS catalogue within 200" 
of the source). 

3. Data analysis 

3.1. Basic parameters 

The observed SEDs have to be characterised in or- 
der to be able to compare them and learn about the 
object properties. We computed Tboi following the 
method described in Chen et al. (1995a). This is 
the temperature of a black body, which posseses the 
same mean frequency as the respective source. The 
bolometric luminosity is determined by integrating the 
SED and adapting a distance. Note that Lboi might 
be underestimated by this method, when the emission 
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maximum is not well covered (no data from 100 to 
350 fim). We compute the ratio of bolometric to sub- 
mm (A > 350 ^m) luminosity. This ratio is used in the 
source classification (see Sec. 3.2). 

The mean spectral index a was determined by a 
powerlaw fit to fluxes measured longward of 400 ^m. 
a is converted to the powerlaw index of the dust emis- 
sivity (assuming optically thin emission) by f3 = a-2. 
The powerlaw fit allows for a uniform determination 
of the flux density at 1.3 mm, even if the object was 
not observed at this wavelength. The fluxes are scaled 
to a distance of 300 pc and transformed into enve- 
lope masses by Mcnv = 1-5 Mq * Fi°3°P;:^[Jy] (follow- 
ing Motte and Andre (2001)). We used a dust opacity 
at 1.3 mm of 0.01 cm^ g~^ and a dust temperature of 
24 K for all objects. A lower dust temperature of 15 K, 
as e.g. suggested by Motte & Andre (2001) for the 
Taurus sources, would lead to envelope masses which 
are a factor of 1.6 higher 

The data for some sources are inhomogeneous. Dif- 
ferent aperture sizes are used at different wavelengths. 
To ensure a consistent treatment of the data for all 
sources, datapoints were selected or excluded manu- 
ally before the determination of the source properties. 
In Table 4 (Appendix B, available online only) we list 
datapoints not used for the determinations, together 
with further notes. The obtained parameters for all our 
sources are listed in Table 1 . This table is divided in 
three sections. First we list all sources where the SED 
is very well sampled and all parameters can be de- 
termined accuratly. The second part contains sources 
where the SED is only well determined on one side 
of the peak and for the other side only upper limits 
(e.g. VLA 1623) or very few data (e.g. CB232) are 
available, hence mostly limits for the source proper- 
ties could be derived. In the last part we list sources 
that have very insufficient data and where at best only 
P and Mcnv or a lower limit for Lboi can be deterimed. 
The objects Trifid-TC 3, GF9-2, and IC 348 MMS are 
not listed in this table since none of the source proper- 
ties could be determined. 

3.2. Source Classification 

To classify the investigated objects we conducted 
three different tests for objects listed in the first two 
parts of Table 1 . A source was negatively tested to be 
a Class object when; (1) Tboi is larger than 80 K; 
(2) Lsnim/Lboi is less than 0.005; (3) There is a NIR 
(A <5 /im) detection. If all three tests are positive then 



Fig. 1. — Distribution of Lboi (filled circles) and Menv 
(open circles) against Tboi of the Class and Class 0/1 
sources with distances of less than 500 pc. 
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the source was marked as Class object (0). If at least 
two tests are positive and the third could not be con- 
ducted or was negative, then the source is marked as 
borderline object (0/1). In case of two negative tests 
the source is classified as Class 1 (1). For all remain- 
ing objects a classification cannot be performed and 
the source is marked by '?'. In principle we further 
need to test if the object possesses an extended en- 
velope. Hence, we marked all sources with a ^ sym- 
bol, where non or no conclusive sub-mm or millimeter 
map is available. The lack of such data mostly effects 
southern hemisphere objects, which anyway could not 
be properly classified by our tests. 

The classification obtained for each source accord- 
ing to these criteria is listed in the last column of Ta- 
ble 1 . Positions in the Tboi-Lboi diagram of all sources 
properly classified as Class or Class 0/1 can be seen 
in Fig. 1 together with the Tboi-Monv positions. In or- 
der to obtain a well defined sample with as few as pos- 
sible selection effects, we excluded sources with dis- 
tances larger than 500 pc in the diagram. 

3.3. Evolutionary Models 

Bolometric temperature, luminosity and envelope 
mass are three observational quantities available for 
these young sources. To investigate physical proper- 
ties, such as final masses and ages, we need to com- 
pare these values with evolutionary models. For the 
comparison we selected all sources from part one in 
Table 1 and the objects from part two with a clear clas- 
sification (0, 0/1, or 1), and where Tboh Lboh and Mcnv 
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could be determined. 

There are several evolutionary schemes available in 
the literature. We use three different models to esti- 
mate ages and final masses of our sources. (1) The 
evolutionary diagram presented in Fig. 12 of Myers et 
al. (1998), where the infall rate matches the isothermal 
sphere infall solution at early times and exponentially 
declines later. (2) Evolutionary tracks shown in Fig. 6b 
in Andre et al. (2000), where accretion rate and enve- 
lope mass decline exponentially. (3) The evolution- 
ary scheme developed by Smith (1998; 2000; 2002), 
where the mass accretion rate first increases exponen- 
tially and then shows a powerlaw fall off. 

All of these evolutionary models are subject of 
shortcommings. They are based on assumptions about 
the time dependence of the mass accretion rates. While 
the overall decline with time has substantial observa- 
tional support (Calvet et al. (2000)), the particular be- 
haviour is unclear Analytical and numerical models 
reproduce both, exponential decline (Shu et al. (2004)) 
or well defined early peaks followed by a fall off 
(Schmeja & Klessen (2004)). The usage of Tboi and 
Lboi as input parameters (Smith (1998; 2000; 2002), 
Myers et al. (1998)) will lead to uncertainties since 
both depend on the viewing angle, in contrast to the 
envelope mass (Andre et al. (2000)), which is inferred 
from optically thin dust emission. Hence, the absolute 
values of the determined ages and masses are model 
dependent but the relative values can be used to put 
these youngest known protostars in an evolutionary 
sequence. Calculated ages and final masses from all 
models are listed in Table 2. 

4. Discussion 

4.1. Source classification 

We investigated SED's of 95 young protostellar 
sources. 59 of these objects (60%) possess sufficient 
observational data to accurately determine all proper- 
ties (TboN LboN Mcnv)- Out of thcsc, 27 objects could 
be unambiguously classified as Class object, 23 as 
Class 0/1, and 9 as Class 1 protostars (according to our 
definition). There is a remaining relatively large num- 
ber (36 objects; 40 %) where we do not have sufficient 
information about the SED for a proper classification. 

There are two main gaps in the observational data 
which mean this large fraction could not be classified 
properly. ( 1 ) Only data at wavelengths shorter than the 
emission maximum and non or only insufficient sub- 



mm or millimeter observations are available. Hence, 
only upper limits for Tboi and lower limits for Lboi can 
be determined. (2) Only data at wavelengths longer 
than the maximum and no observations or just upper 
limits at shorter wavelengths are available. For these 
objects only f3 and Mcnv can be estimated properly. 

4.2. Uncertainties 

All determined object parameters are subject to dif- 
ferent errors, due to their calculation. The powerlaw 
index of the opacity is determined from the mean slope 
of the SED at A > 400 fim. Using fluxes taken with dif- 
ferent aperture sizes to determine (3 will lead to differ- 
ent results. Estimating /3 as the mean slope ensures that 
for each object a 'mean' aperture size for the (sub)mm 
data is used. Typical errors of the estimated values are 
0.3. 

The luminosity is determined by integrating the 
SED and adopting a distance. An error of 20% for the 
distance (a typically value for most of the objects with 
d < 500 pc) leads to 40% uncertainty for the luminos- 
ity. Also the luminosity might be underestimated when 
the SED is not well sampled at the emission peak, or 
when huge parts of the SED are only described by up- 
per Umits (e.g. VLA 1623 orHH24MMS). Additional 
uncertainties are expected in cluster regions, where 
current far-infrared instruments do not possess suffi- 
cient angular resolution to separate individual sources 
(e.g. SVS 13 B). Non-detections or observations in the 
NIR or at millimeter wavelength have almost no influ- 
ence, since most energy is radiated near the emission 
maximum. All together the determined source lumi- 
nosities are uncertain by about 50% in most cases. 

For the bolometric temperature, datapoints at short 
wavelengths are important. In the case of well sampled 
SEDs from the NIR to the sub-mm, Tboi can be deter- 
mined with about 5 to 10 K accuracy. In the case of 
upper limits in the NIR, the same accuracy applies as 
long as the limits are at least 5-6 orders of magnitude 
lower than the flux at the maximum of the emission. 
Note that the sources classified as Class 0/1 have typ- 
ically NIR detections 4-5 orders of magnitude below 
the emission maximum, while for the Class sources 
we find typical non-detections 5-7 orders of magnitude 
below the maximum. In case of 'bad' NIR limits the 
uncertainty in Tboi can well exceed lOK. Hence for 
our Class sources the temperatures are accurate to 
10 K, while for the Class 0/1 objects they are certain to 
20 K. Note that for older sources (Class 0/1 or Class 1) 
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viewing angle effects might become important. They 
lead in case of edge-on sources to underestimates of 
Tboi and Lboi- Subsequently the Lsmm/Lboi ratio is 
overestimated, influencing the source classification. 

The envelope mass is also subject of various 
sources of errors. Beside the problems with measuring 
the fluxes in different apertures, the distance is used in 
the determination. Further the applied dust opacities 
are uncertain by a factor of two (e.g. Motte & Andre 
(2001)), and the dust temperatures are not very well 
known. Hence, the masses might suffer from errors of 
a factor of three. 

How do these errors influence inferred ages and 
masses? The models of Smith (1998; 2000; 2002) and 
Myers et al. (1998) use Tboi and Lboi as main input 
parameter In both cases the error in the measured 
luminosity leads to uncertainties of about a factor of 
two for the estimated final masses. Errors in Tboi will 
transform into age uncertainties which are about 30%. 
The model of Andre et al. (2000) uses the envelope 
masses as input. Hence inferred final masses are un- 
certain by a factor of three and estimated source ages 
are uncertain by a factor of two. Note that these un- 
certainties do not take into account systematic errors 
due to the various assumptions made in the different 
evolutionary models. 

4.3. Distribution of source parameters 

To investigate the statistical properties of the ob- 
tained Class sample, we restrict all the following dis- 
cussions to sources classified as Class or Class 0/1 
objects within 500 pc. 

The power-law index of the opacity is evenly dis- 
tributed from 0.0 to 2.0 with a broad maximum be- 
tween 0.5 and 1.5. For our sources Tboi ranges from 
30 to 90 K, where the 90 K limit is due to our classi- 
fication criterion for Class sources. The number of 
sources gradually rises from 30 to 70 K and falls off at 
higher temperatures. One expects a rising (or at least 
constant) number of sources towards higher temper- 
atures, because with higher temperature an increased 
age and a slower evolution in temperature is predicted 
by the models. Hence, the Tbd distribution shows that 
there are missing objects in our sample in this tem- 
perature range. This might be due to a classification 
of these objects as Class 1 protostars and reflect the 
fact that sources with the same bolometric tempera- 
tures are not necessarily in the same evolutionary state. 
The distribution of Lbd corresponds very well with the 



mass function (see below), due to the assumptions in 
the evolutionary models. Similarly the distribution of 
Monv shows a peak at about one solar mass. At lower 
masses our sample certainly suffers from incomplete- 
ness. The Lgnim/Lboi ratio shows a broad distribution 
between 0.01 and 0.1. Only very few larger values are 
found (IRAM 04191, L 1448 NW). 

We further investigated if the source properties are 
correlated. There is no correlation between the dust 
opacity f3 and the parameters Tbd and Lboi (correla- 
tion coefficients (c.c.) -0.034 and -0.102). A weak 
correlation (c.c. -0.412) is found between f3 and Mcnv, 
but this is due to a few high envelope mass sources. 
Also Tboi is not correlated with Lbd and Monv (c.c. 
0.149 and -0.219). On the other hand the envelope 
mass shows a correlation with the bolometric luminos- 
ity (c.c. 0.660). This correlation indicates that objects 
with more massive envelopes possess a higher lumi- 
nosity and hence higher accretion rates. 

4.4. The different evolutionary models 

Due to the different model assumptions about the 
time evolution of the mass accretion rates, all three in- 
vestigated evolutionary models lead to different values 
for ages and final star masses. The ages especially are 
very model dependent. We investigate if there are cor- 
relations between the inferred values for age and fi- 
nal mass for the different evolutionary models in our 
Class and Class 0/1 sample. 

On first sight the ages differ significantly between 
the models. Ages from Myers et al. (1998) are much 
larger than usually assumed for Class sources (a cou- 
ple of 10^ yrs). In the model of Andre et al. (2000) the 
evolution in the first 10000 yrs is very fast, and hence 
the inferred ages are on average very smafl (< lO** yrs). 
However, there is a weak correlation between the ages 
inferred from the model of Smith (1998; 2000; 2002) 
and Myers et al. (1998) (c.c. 0.425). Comparing the 
models from Andre et al. (2000) with Myers et al. 
(1998) and Smith (1998; 2000; 2002) we find no cor- 
relation (c.c. 0.143,0.263). 

The predicted final masses are much more impor- 
tant than the ages, since the mass distribution should 
represent the observed initial mass function. Here we 
investigate if the relative masses obtained from the dif- 
ferent models are comparable. There are no obvious 
differences in the inferred masses between the mod- 
els. The models of Myers et al. (1998) and Andre et 
al. (2000) are limited in their mass range (0.3-0.7 and 
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0.2-3.0 Mq, respectively). We find good correlations 
between masses obtained from the different models: 
c.c.=0.676 for Smith (1998; 2000; 2002) and Myers et 
al. (1998); c.c.=0.656 for Smith (1998; 2000; 2002) 
and Andre et al. (2000) and c.c.=0.771 for Myers et al. 
(1998) and Andre et al. (2000)). 

4.5. The Mass Function 

The evolutionary scheme by Smith (1998; 2000; 
2002) allows the comparison of all three main obser- 
vational quantities (Tboh Lboi, Menv) with the model. 
Also there is no restriction on the final mass as in 
the other models. Further, the assumed mass accre- 
tion rates are very similar to accretion rates obtained 
by hydrodynamical simulations of star formation (e.g. 
Klessen (2001), Schmeja & Klessen (2004)). Hence it 
is worth to investigate the resulting final mass function 
obtained for our sources using this model. Certainly 
the absolute values for the masses are uncertain, but 
giving the good correlations to the other models (see 
above), the relative masses, needed to determine the 
slope in the mass function, are correct. 

In Fig. 2 we present the resulting mass func- 
tion for our sample of Class and Class 0/1 objects 
within 500 pc. The slope determined for objects with 
M > 0.5 M0 is -0.9 ± 0.2. Considering the low num- 
ber of sources, and that no binary correction is done, 
this is in good agreement with the Salpeter slope of 
-1.35 measured for the solar neighborhood. We esti- 
mate a completeness limit in our sample for objects 
with final masses below 0.4 M0. 

There is, however, a number of predicted very 
low mass sources. These are objects where the 
model of Smith (1998; 2000; 2002) seems not to 
work properly. Looking in detail, we find that 
the evolutionary scheme of Smith (1998; 2000; 
2002) can explain the all observed properties (Tboi, 
Lboi, Menv) for most of the sources simultaneously. 
In case of eleven objects, however, the scheme is 
not working (IRAS 03256H-3055, IRAS 03282H-3035, 
B213, IRAM 04191, L1448N, IRAS 04325H-2402, 
IRAS 04368H-2557, IRAS 15398-3359, HH24MMS, 
NGC 2068 LBS 17, VLA 1623). These sources pos- 
sess a much lower luminosity (considering their Tboi 
and Menv values) than the average of the objects. 
There are also three objects where the luminosity is 
much above the average (NGC 1333-12, HH212-MM, 
L 1641 N). Note that this applies also for the more 
distant objects L 1246-SMM 1 (luminosity lower) and 



Fig. 2. — Mass Function for the properly classified 
Class and Class 0/1 objects within 500pc, obtained 
from the model of Smith (1998; 2000; 2002). 
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CepE, IRAS 20050H-2720, and MonOBl IRAS12 SI 
(luminosity higher). For the remaining sources the 
model explains very well all three observed proper- 
ties. 

There are two reasons for this inabiUty to explain 
the data of all our objects. (1) One might be found in 
the mass accretion rates obtained by Klessen (2001) 
and Schmeja & Klessen (2004). It turnes out that 
their accretion rates on average show a time evolu- 
tion similar to the one used in the model of Smith 
(1998; 2000; 2002). But on short timescales the accre- 
tion rate varies significantly from this average. There 
are periods with lower and higher accretion rates re- 
sulting in higher and lower luminosities. Hence, the 
very low luminosity sources might represent a quiet 
accretion phase, while the higher luminosity sources 
represent phases of enhanced accretion. (2) A second 
reason might be that these sources do not accrete their 
mass in the same way. Hence, the mass accretion rate 
shows a signifficantly different evolution in these ob- 
jects. 

In case of the latter, this suggests two different 
groups of sources. 'Normal' sources, possessing ac- 
cretion rates similar to that obtained by Schmeja & 
Klessen (2004), and 'abnormal' sources showing on 
average lower luminosities/accretion rates, consider- 
ing their Tboi and Monv values. Note that viewing 
angle effects could be responsible for some of these 
objects. These low luminosity sources might repre- 
sent objects where ambipolar diffusion dominates the 
accretion process instead of turbulence. A detailed 
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analysis and discussion of this subject can be found 
in Froebrich et al. (2004). 

4.6. The different Star Forming Regions 

Does the star formation process, in particular the 
time dependence of the mass accretion rate vary with 
the star forming region or is it a uniform function? 
Due to the low number of sources a proper statistical 
analysis of this question is not possible yet. However, 
in order to see possible trends we chose four regions 
(Perseus, Taurus, Orion, and Serpens) where we have 
a sufficient number of confirmed Class or Class 0/1 
objects (13, 4, 10, 6, respectively) and compared these 
regions with the average of all sources. 

In our sample of Class and Class 0/1 sources about 
25% do not follow the Tboi-Lboi-Mcnv relation of the 
majority of our objects. These sources possess much 
lower luminosities than suggested by their bolometric 
temperature and envelope mass. Hence, they might be 
in a phase of lower mass accretion or are gaining their 
mass in a different way. Perseus reflects the overall 
average with three out of 13 objects. The same ap- 
plies for Orion (two out of ten), where we might miss 
low luminosity objects due to the larger distance. In 
Serpens all of the six sources show a 'normal' be- 
haviour. This is, according to the small number of 
objects, still in agreement with the average. Taurus, 
however, shows a very different picture. Here all four 
Class or Class 0/1 objects fall into the low luminos- 
ity category. Objects showing a 'normal' behaviour 
should have been easily detected since their higher lu- 
minosity. Hence, even if there are just four sources this 
strongly suggests that in a region of isolated star for- 
mation (as in Taurus) the time evolution of the mass 
accretion rates is different from regions where stars 
form in clusters (in agreement with previous findings 
from Henriksen et al. (1997)). Due to the bad statis- 
tics nothing can be said about differences among the 
clusters. 

5. Conclusions 

A literature search for photometric broadband 
observations of 95 confirmed or candidate Class 
sources was conducted. To our best knowledge all 
available broadband photometric data was used to con- 
struct SEDs from I fim to 3.5 mm. If possible we 
determined basic properties of the sources (sub-mm 
slope of the SED, Tboi, Lboi, Lsmm/Lboi, and Mc„v). 
For 59 objects sufficient enough data are available for 



a proper determination of the source parameters. 27 of 
these are classified as Class 0, 23 as Class 0/1, and 9 as 
Class 1 protostars. 

To investigate the statistical properties of the ob- 
tained sample of very young protostars we used all 
objetcs within 500 pc which are properly classified as 
Class or Class 0/1 to determine age and final star 
mass. Therefor we used three different evolutionary 
models for protostars from Smith (1998; 2000; 2002), 
Myers etal. (1998), and Andre et al. (2000). Consider- 
ing the uncertainties in the measured source properties 
the predicted final star masses from the different mod- 
els show a good correlation. The absolute ages are, 
however, very model dependent. 

An investigation of the final star masses and the re- 
sulting mass function shows a good agreement with the 
IMF for stars with M > 0.5 Mq. Our obtained Class 
sample is limited to sources with final masses above 
0.4 Mq . A number of objects (25 %) is found that pos- 
sess a much lower luminosity than the rest of our sam- 
ple, considering their bolometric temperature and en- 
velope mass. These objects might be in a 'quiet' accre- 
tion phase or the time evolution of their mass accretion 
rate is significantly different from the majority. This 
might especially be the case in the Taurus star forming 
region, where all identified Class or Class 0/1 sources 
belong to this group. 
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Table 1:: Obtained object parameters, sorted by SED quality 
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Table 1 - continued from previous page 
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1.2 








0.49 






G 34.24+0. 13 MM** 


0.9 








2.8 






RCrARl** 






>3.5 










IRAS 19180+1114** 






>5.2 










S 106-SMM** 


1.1 








4.9 






IC1396W** 






>26 










NGC 7538 S** 


2.2 








51 







° Object name used in this paper. Further common names are given in Table 3. The * markes objects where the SED 
is is only well determined on one side of the peak and for the other side only upper limits or very few data is available. 
The ** markes objects with too few observational data to properly estimate even a limit for Tboi- 
^ Classification due to the procedure described in Sect. 3.2. In brackets we give as 'y', 'n', or '?' the answers for our 
three different tests if the object is a Class source. The ending number gives the orders of magnitude between the 
maximum of the emission and the NIR detection/upper limit. 

t Classification uncertain, since non or no conclusive sub-mm or millimeter maps are available. 
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Table 2:: Age and Mass estimates from evolutionary models 



Object" 




Class'' 




;e'=[I03 yrs] 




M|„ [Mq] 










Smith 


Myers 


Andre 


Smith 


Myers 


Andre 


L 1448 NW 





(yyy6) 


23 


/ J 


c 

J 


0.2 


0.8 


1.2 


L1448C 





(yyy6) 


37 




o 
o 


0.4 


1 


1.1 


RNO 15 





(yyy6) 


47 


LI J 


1 AM 
iOU 


0.6 


>2 


2.0 


IRAS 03256+3055 


ot 


(yyy5) 


36 




^M 


0.05 


0.4 


0.4 


NGC 1333-12 





(yyy7) 


33 


^MM 


A 

4 


2 


>2 


>3 


NGC 1333-14 A 





(yyy6) 


29 


OMM 





1 


>2 


>3 


NGC 1 333-14 B 





(yyy6) 


29 


OMM 





1 


>2 


3.0 


IRAS 03282+3035 





(yyy5) 


37 


1 1 M 


o 
O 


0.06 


0.4 


0.7 


HH211-MM 





(yyy6) 


25 


/ J 


Q 
O 


0.3 


0.8 


0.8 


IRAM 04191 





(yyy5) 


24 


9M 


M 


0.01 


<0.3 


0.5 


OMC3-MM6* 





(yyy7) 


48 




o 
o 


5 


>2 


>3 


HH212-MM 





(yyy5) 


35 


i /j 


AM 


0.4 


1 


0.45 


HH25 MMS 





(yyy5) 


28 


1 ^M 


-7 
/ 


0.4 


1 


1.3 


HH24MMS* 





(yyy5) 


53 


TMM 





0.8 


>2 


3.0 


NGC 2068 LBS 17* 





(yyy5) 


27 


/ 3 


M 


0.10 


0.4 


2.5 


NGC 2264G-VLA2 





(yyy5) 


44 


250 


25 


0.6 


>2 


0.7 


VLA 1623* 





(yyy7) 


36 


150 


6 


0.1 


0.5 


0.9 


IRAS 16293-2422 





(yyy7) 


29 


200 


5 


1 


>2 


>3 


IRAS 16544-1604 





(yyy5) 


39 


150 


10 


O.Oo 


A /I 

0.4 


A 1 

U.3 


Serp-S6s N 


u 


(yyy5) 


5j 


150 


6 


U.z 


A O 


i.i 


Serp-i)MM4 


u 


(yyy6) 


o o 
3o 


200 


6 


U.z 


A O 


i.J 


Serp-SMM i 


u 


(yyy5) 




200 


8 


U.z 


A O 

U.o 


1 A 
i.U 


Serp-SMM / 





(yyy5) 


47 


200 


40 


A 'I 

0.2 


A T 

0.7 


A /I 

0.4 


IRAS 19345+0727 


r\ 




(yyy6) 


29 


125 


8 


A O 

O.z 


A *7 

0.7 


A O 

0.8 


IRAS 203o6+6/51 


u 


(yyy5) 


z9 


150 


8 


U.J 


i.i 


1 o 
i.z 


L, 944-SMM I 


f\ 


(yyy5) 


C 1 


200 


18 


A O 1 

U.zi 


A T 


A C 
U.J 


Cep b-MM 


u 


(yyy8) 


3o 


400 


20 


4 


>z 


z.o 


IRAS 23238+7401 





(yyy5) 


36 


100 


20 


A A /I 

0.04 


A A 

0.4 


^A T 

<0.2 


L 1448-12 


0/1 


(yyn6) 


29 


1 ^M 


O 

o 


0.4 


0.9 


1.1 


L 1448 N 


O/I 


(yyn4) 


44 







0.4 


1 


2.1 


RNO 15 FIR 


0/1 


(yyn5) 


39 


225 


30 


0.4 


1 


0.7 


NGC 1333 II 


0/1 


(?yy6) 


74 


325 


60 


1 


>2 


0.8 


B2I3 


0/1 


(yyn4) 


41 


90 


5 


0.01 


<0.3 


0.2 


IRAS 04325+2402 


0/1 


(yyn4) 


43 


150 


30 


0.04 


0.4 


0.7 


LI527 


0/1 


(yyn4) 


34 


125 


5 


0.1 


0.6 


0.9 


HHI14MMS* 


0/1 


(?yy6) 


70 


350 


10 


1 


>2 


1.6 


IRAS 05173-0555 


0/1 


(?yy5) 


130 


250 


30 


0.7 


0.8 


0.5 


RN0 43-MM 


0/1 


(yyn5) 


35 


200 


20 


0.3 


0.9 


0.6 


LI64IN 


O/I 


(yyn5) 


43 


>500 


9 


4 


>2 


>3 


HH 1 1 1 MMS 


0/1 


(yyn5) 


54 


300 


20 


1 


>2 


1.2 


MonOBlIRAS12Sl 


0/1 


(yyn5) 


20 


150 


8 


7.8 


1.3 


>3 


IRAS 15398-3359 


O/I 


(yyn4) 


36 


100 


7 


0.04 


0.4 


0.4 


IRAS 18148-0440 


0/1 


(yyn5) 


34 


200 


20 


0.5 


1 


0.8 



Continued on next page 
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Table 2 - continued from previous page 



Obiect" 




Class'' 




;e^[103 yrs] 














Smith 


Myers 


Andre 


Smith 






Serp-FIRS 1 


0/1 


(yyn6) 


32 


300 


8 


2 


>2 


>3 


IRAS 18331-0035 


0/1 


(yyn4) 


31 


150 


9 


0.2 


0.7 


0.8 


IRAS 20050+2720 


0/1 


(yyn5) 


53 


>500 


10 


10 


>2 


>3 


IRAS 2 1017+6742 


0/1 


(yyn5) 


43 


125 


40 


0.04 


0.4 


<0.2 


CB230 


0/1 


(yyn4) 


43 


225 


17 


0.34 


0.9 


0.6 


NGC 7129 FIRS 2* 


0/1 


(yyn4) 


49 


>500 


8 


14 


>2 


>3 


L1246-SMM1 


0/] 


(?yy5) 


290 


275 


8 


1.93 


0.8 


1.0 


IRAS 23385+6053* 


0/] 


(yny7) 


36 


>500 





>50 


>2 


>3 


L1551-IRS5 




(nyn4) 


290 


350 


9 


6 


>2 


1.6 


L1551-NE 




(nyn5) 


250 


225 


10 


1 


0.7 


0.7 


OMC3-MM9 




(nnn6) 


370 


>500 


20 


50 


>2 


2.8 


OMC3-MM7 




(nyn4) 


380 


>500 


30 


20 


>2 


1.7 


HH147MMS 




(nnn3) 


430 


>500 


150 


30 


>2 


2.0 


IRAS 08076-3556* 




(n?n5) 


365 


300 


200 


4.2 


1.0 


2.5 


SerpSMM5 




(nyn4) 


350 


250 


5 


1 


0.6 


1.5 


IRAS 18148-0440 




(nyn3) 


280 


250 


8 


0.9 


0.6 


0.8 


L1251B 


l1 


(nyn4) 


61 


250 


7 


0.4 


1.0 


2.0 


L1211* 


1 (nyn4) 


426 


>500 


8 


33 


>2 


1.2 



" Object name used in this paper. Further common names are given in Table 3. The * markes objects where the SED 
is is only well determined on one side of the peak and for the other side only upper limits or very few data is available. 
^ Classification due to the procedure described in Sect. 3.2. In brackets we give as 'y', 'n', or '?' the answers for our 
three different tests if the object is a Class source. The ending number gives the orders of magnitude between the 
maximum of the emission and the NIR detection/upper limit. 

Ages and final masses obtained from the evolutionary models of Smith (1998; 2000; 2002), Myers et al. (1998), and 
Andre et al. (2000). 

^ Classification uncertain, since no or no conclusive sub-mm or millimeter maps are available. 
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Source Sample 



Table 3:: Complete source sample 



Object" 


Ref.Obj.'' 


RA(2000) 


Dec (2000) 


Ref.'^ 


d[pc] 


Ref.'^ 


W3 OH 


15,43 


A'^ T7 A/1 

Oz 21 04.72 


1 /c 1 A n 

+ol 52 24.7 


71 


OO AA 

2200 


A T 

43 


TW3(H20) 














L 1448-12 


22, 43 


AT C 

03 22.D 


1 T A /I c r\/z 

+30 45 Oo 


C A 

59 


T A A 

300 


C A 

54 


IRAS 03222+3034 














L 1448 NW 


48, 54, 70 


AT C T C /C C 

03 25 35.65 


1 T A /I C T /I 

+30 45 34.2 


A n 

47 


T A A 

300 


C A 

54 


LDN 1448 IRS 3C 














L 1448 N 


16, 22, 40, 43, 48, 54 


AT C T /C T 

03 25 36.3 


1 T A AC 1 /C 

+30 45 16 


A O 

48 


T A A 

300 


C A 

54 


LDN 1448 IRS 3, IRAS 03225+3034 












T 'X A AO 


lo, 22, 40, 43, 4o, j4, /U 


AT O^T "2 C Q 

U3 zj 3o.o 


+30 44 00 


A 

4o 


300 


<A 

54 


L 1448-MM 














KJNU i J rlK 


AO £,A -TA 

46, 04, /U 


AT TA A 

03 27 39.0 


1 TA 1 O CA 

+30 12 59 


*7 1 

71 


T C A 

350 


*7A 

70 


LDN 1455 FIR 














KNO 15 


TA /n\ 

70 (?) 


AT on /1 Q 
Oj z / 4z.y 


+30 ill! 


IK) 


3jO 


IK) 


L 1455, B 204 














T<lf~^i~^ 1 111 T1 

JNuC 1333 11 


70 (ClassO/1) 


AT OQ QO 1 

03 zo 3o. / 


1 T 1 1 T TO 

+31 13 3Z 


n 1 

/ 1 


350 


^A 

54 


IRAS 03255+3103 














IRAS 03256+3055 


68 (Class 1) 


AT A A 

03 Z5 44.5 


1 T 1 AC TA 1 

+31 05 39. / 


/; o 
00 


A AA 

400 


2 


XT/""/"' 1 TOT T'^ 

JNuC 1333-12 


40, 43, 54, /O 


AT OQ CC /I 


1 T 1 1 /l T C 

+31 14 3j 


JO 


3jO 


KA 

j4 


IRAS 03258+3104 














SVS 13 B 


//^r^x /IT /I O i \ /I A 

22 (OF), 43, 48 (Classl), 49 


AT '^A AT A/; 

03 zy 03.00 


1 T 1 1 C C 1 T 

+31 15 51. / 


A n 
4/ 


T CA 

350 


TA 
/O 


INOC 13jj-I4A 


Ar\ A1 <A 

2.2., 4U, 4j, j4 


AT OQ 1 r\/i 
03 Zy iZ.04 


1 T 1 IT T A C 

+31 13 30. J 


IK) 


3jO 


j4 


XT/"'/"' 1 TIT T/1 1 ) 

INCjC 1333-I4I3 


22, 43, 54 


AT '^A IT/; 

03 zy 13.0 


1 T 1 IT A/; /; 
+31 13 Oo.o 


Ik) 


T CA 

350 


C A 

54 


IRAS 03282+3035 


l/C /I") AO C/1 

16, 43, 48, 54 


AT T 1 OA T 

03 31 Z0.3 


1 T A /I C O C 

+30 45 25 


A O 

46 


T A A 

300 


C /I 

54 


IC 34o MMi) 


63 


AT /I T c/r A 

03 43 JO.V 


1 "JO r\'^ c\f^ 
+32 03 Oo 


o3 


300 


<A 

54 


T TT T'^11 \ /f\ /f 

HH 211 -MM 


/f T C A HA '~l f\ 

43, 54, 64, 70 


AT /IT C/C 

03 43 56 


1 T O A A /I O 

+32 00 48 


C A 

54 


T AA 

300 


C A 

54 


B 213 


16 (OF), 48 (Classl), 


r\A 1 A /I T A A 

04 19 43.00 


1 0*7 1 T T T T 

+27 13 33.7 


/I A 

60 


'i A n 

140 


C A 

54 


IRAS 04166+2706 


54 (O/I), 68 












IR AM 04191 


43,54 


A/1 1 C/T A 

04 Zl JO. 9 


1 1 ^ oo A A 

+ 15 29 4o 


3o 


1 /I r\ 


^A 

54 


L1551-1RS5 


13, 22 (OF) 


A/1T1 T/1 IC 

04 31 34.15 


1 1 O AO AC O 

+ 18 08 05.2 


A A 

44 


140 


C A 

54 


IRAS 04287+1801 














L1551-NE 


13,22 


f\A T 1 /I /I /I /I 

04 3 1 44.44 


1 1 r\c n 
+ io Oo J/ 


3 / 


1 /I r\ 
14U 


<A 

j4 


IRAS 04325+2402 


26 (Tboi = 157K) 


A /I T C T C A 

04 35 35.0 


1 O/I AO O O 

+24 08 22 


54 


140 


2 


LDN 1535 IRS 














L1527 


13, 16, 43,48,49, 54 


A/1 TA CT A 

04 39 53.9 


1 O/C AT 1 1 

+26 03 1 1 


A O 

48 


1 /I A 

140 


C A 

54 


IRAS 04368+2557 














HH114MMS 


24, 43 


AC 1 O 1 C 1 

05 18 15.21 


1 A*7 1 O AT 

+07 12 03 


o 1 
21 


A CA 

450 


A T 

43 


IRAS 05173-0555 


70 


05 19 48.9 


-05 52 05 


61 


460 


70 


RNO40FIR, L1634 














LI 634 IRS 7 


70 (ClassO?) 


05 19 51.5 


-05 52 06 


70 


460 


70 


RN0 43-MM 


7, 43 


05 32 19.4 


+ 12 49 41 


32 


400 


43 


IRAS 05295+1247 














OMC3-MM6 


25,43 


05 35 23.4 


-05 01 29 


41 


450 


43 


OMC3-MM9 


25 


05 35 25.9 


-05 05 44 


41 


450 


43 


OMC3-MM7 


25 


05 35 26.5 


-05 03 55 


41 


450 


43 


IRAS 05329-0505 














HH1-2MMS3 


49 


05 36 18.2 


-06 45 45.3 


49 


460 


49 


L1641N 


22 (OF), 70 


05 36 18.6 


-06 22 10 


70 


390 


70 


IRAS 05338-0624 














HH1-2MMS2 


49 


05 36 18.8 


-06 45 25.3 


49 


460 


49 


L1641-VLA1 


43,49 


05 36 22.8 


-06 46 07.6 


49 


460 


49 


HH1-2MMS 1 














L 1641 SMS 111 


70 (ClassO/ 1) 


05 36 24.0 


-06 24 54 


70 


390 


70 


HH147MMS 


49 


05 36 25.2 


-06 44 39.8 


49 


460 


49 



Continued on next page 
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Table 3 - continued from previous page 



Objecf" 


Ref.Obj.'' 


RA (2000) 


Dec (2000) 


Ref.'^ 


d[pc] 


Ref. 
















T 1 f\A^ 'sM's TV 




OS "^^i A^ 8 
UJ J\J +1.0 


OA 9A 1 9 

— UO ZO IZ 


70 

/ u 




70 

/ u 


MnP 709'^ MM 1 


1Q A"^ 


ns A^ 7A szL 


09 1 8 00 

— UZ lo uy 


so 

jy 


Asn 

4jU 


AT. 
4J 


INvJv^ ZUZ^H—rUV J 


U, 4j 


OS /L1 /L/L 99 
uj> +1 ++.ZZ 


01 Al ^9 
— Ul JJ 41 . JZ 


AS 

OJ 


Asn 

4JU 


A'K 
4J 


INO^^ ZUZ4-rlK0 


A /I 




01 ^A 00 ^A 
— Ul JO UU.jO 


AS 
OJ 


4jU 


A 7 
4J 


1414 9 1 9 MM 


7 A"^ 
/, 4 J 




01 0"^ 01 
— Ul UJ Ul 


71 

/ 1 


Ann 

4UU 


A'? 
4J 


IKAo Uj41 J-U1U4 














1414 9^ MMQ 

ori z J iviivio 


1 9 A'^ 
iZ, 4 J 


AC /lA 07 


00 1 A 1 
— uu IJ 41 


JZ 


A^n 

4JU 


A'K 
4J 


TJLJ '7/1 \yr\yrc 
rlo z4 iVIiVIo 


99 Al. 
zz, 4j 


/lA MQ Q 


00 1 /1 7 
— UU lU 4 / 


S9 

JZ 


A 

4jU 


A 'K 
4J 


INOV^ ZWOo L/Do i / 


4j 


A^i '?0 R 
UJ 40 jV.o 


00 09 AO 
— UU UZ 4U 


SA 

JO 


Ann 

4UU 


AS 
4J 


TJLJ 1 1 1 \yr\yrc 
orl ill iVllVlo 


4V 


Uj J 1 40. J 


_l09 /I 8 98 
+Uz 40 ZO 


AQ 

4y 


/I An 

40U 


AQ 

4y 


IKAo U j4y 1 +UZ4 / , 


LJJIN lol / 1 












Mr\r»OU 1 TP A Q 1 9 Q 1 


AQ 


AA Ai R 

UO 41 Uj.o 


1 AQ OA AQ 

-huy J4 uy 


AO 

oy 


snn 


AO 

oy 


T\jnr'99AAn \/t a 9 

IN KJK^ ZZ04 v J- V Z 


14, 4j 


OA A 1 1 Q 

UO 41 lu.y 


1 AQ CA 09 

-huy JO UZ 


1 A 
14 


snn 


A'K 
4J 


IKAo UDJ04+U7JO 
















1 1 A'K 
1 1 , 4 J 


AO AQ TO 
UO UV JZ.O 


TA 0^ 00 

— JO UJ uu 


7 1 
/ 1 


Ann 

4UW 


A'K 
4J 


DTTD 1 
DrlK IZ 














PpH 1 1 nTT?*^ 1 n 

l_-CU i i W livo i \J 




I 1 Ofi ^9 

I I uu JZ 


77 9^ A9 

/ / ZJ T-Z 


90 
zu 


1 sn 


5^ 

J J 


cna-iviivio 1 














C^f^A 1 1 n TP ^ zL 


9A 
zO, 


1 1 flA A7 1 A 
i i UO 4 / . 14 


77 99 "XA 
— / / ZZ J4.U 


JJ 


1 ^n 

1 JU 


S'^ 
J J 


TD A C 1 1 1 H(\(^ 














D n i\ / i-iviivi 


9"^ A'^ 

Z^, 4J 


19m T.1 
IZ Ul J 1 


A^ 08 SA 

— OJ UO J4 


7 1 
/ 1 


9nn 

ZUW 


A'K 
4J 


TD A C 1 1 ^On A/1 ^0 














TP A Q 1 9^nn 7A^s 


9A 
ZO 


1 9 CO TO Q 

IZ jj jo.y 


77 1 ^ ^'^ 

— / / 1 J JJ 


S7 

J / 


1 sn 

loU 


A 



TP A^ 1 9S^^ 7f^^9 


9^1 

zu 


1 Z J / WU. 1 


If^ A8 ^5 

/ 40 J J 


59 

JZ 


1 sn 

1 oU 


A 



TP AC 1 9^'^'i 7A'^ 1 


9A 
ZD 


1 9 0^ ^ 

iz jy uj.j 


77 07 

— 1 1 \} 1 J4 


71 
/ 1 


1 Qn 

loU 


^0 
JU 


TP A C 1 9'^'^/l 7A^^ 


9A 
ZO 


1 9 1 ^ /I 

iz jy 1 J.4 


7A ^1 11 

— /O J 1 11 


71 
/ 1 


1 en 

loU 


A 



TP A Q 1 ^n^A lf,AA 


9A A'K 
ZO, 4j 


1 07 "^A 
1 J U/ JO.U 


77 00 OA A 
— / / UU U4.4 


90 

zy 


onn 

ZUW 


A'K 
4J 


DrlK oO 
















9A A8 
zO, 4o 


1 ^ A'^ 01 
1 J 4 J Ul . J 


-!A AQ 1 9 

— J4 uy iz 


90 

zy 


1 ^n 


A8 
4o 


D ZZo 














TP A C 1 AMI 7 '^Q'^A 


9A 
zO 


1 A 0^ O/l 7 
10 Uj U4. / 


1Q /I ^ O/l 
— jy 4j U4 


7 1 
/ 1 


1 '^n 
1 


9A 
ZO 


\/T A 1 A9^ 

V iOZ J 


1 A 99 A'^ 
10, ZZ, 4^ 


1 A 9A 9A '^9 
10 ZO ZO.JZ 


9A 9A "^0 1 
— Z4 Z4 JU. 1 


A7 
4 / 


1 An 

lOU 


A'K 
4J 


TPA^; 1A9Q'^ 9zL99 
lK/\o iOZ^j-Z^ZZ 


1 a 1 A 99 A'^ 
1 J, 10, ZZ, 4J 


1 A "^9 99 8 
10 JZ Zz.o 


9A 98 T.'^ 
— Z4 Zo JJ 


7 1 
/ 1 


1 An 

lOU 


A'K 
4J 


T 1 ASQ 














TP A C 1 f\AQ'K AOAO 


9A 
ZO 


1 A S9 SA A 
10 JZ JO.O 


A9 A7 SA 

— 4Z 4 / JO 


7 1 
/ 1 


1 J\J 


9A 
ZO 


TP A^ 1 f^SzLd 1 f,OA 




1 fi 57 1 Q 5 

lu J / 1".J 


1 00 9 1 
10 uy zi 


^0 

jy 


1 f^n 
low 


98 
ZO 


T 1 /I A r^T3 AQ 
L 140, V^rS Do 














TrifiH TP ^ 

irinu- 1 J 


1 A'K 
J 1, 4 J 


1 8 09 07 
1 UZ U / 


90 A« 1 1 

— ZJ UJ 1 1 


7 1 
/ 1 


1 Asn 

lOoU 


A'K 
4J 


TP A <i 1 Si 1 ASi C\AAa 
lK/\o ioi4o-W44W 


AH A'^ A8 '^A A9 
4W, 4j, 4o, J4, OZ 


18 17 90 8 
lol / zy.o 


OA '?0 '?8 
— U4 jy Jo 


71 
/ 1 


9nn 
zuu 


SA 
J4 


T AQ^ 














ocip-ouo i> 


1 8 AO A^ 

10, H"W, H- J 


1 8 7Q A8 1 
1 Z^ HO . 1 


-1-01 IfiAl 
\\J 1 1 U H 1 


^8 
JO 


^1 n 


A'K 

4J 


C^rn PTPQ 1 


1 /in A'K 

lo, 4W, 4j 


1 8 90 AO 70 

lo zy 4y. /y 


-1-01 1 ^ 1 8 A 
-hUl IJ lo.O 


'?8 
Jo 




A'K 
4J 


TRAS 1S?7^+01 n 














SerpSMMS 


40 


182951.1 


+01 16 36 


38 


310 


43 


Serp-SMM4 


18,40, 43 


18 29 56.5 


+01 13 10 


38 


310 


43 


Serp-SMM3 


18,40, 43 


18 29 59.2 


+01 13 58 


38 


310 


43 


Serp-SMM2 


18,40 


18 30 00.2 


+01 12 57 


38 


310 


43 


IRAS 18331-0035 


40, 49 


18 35 42.0 


-00 33 18 


71 


310 


24 


HH 108/109 IRS, L 


558 












HH108 MMS 


49 


18 35 46.55 


-00 32 41.8 


71 


310 


24 


G 34.24+0. 13 MM 


34, 43 


18 53 21.4 


+01 13 45.3 


34 


3700 


43 


RCrARl 


22 


19 01 55.8 


-36 57 29.6 


22 


130 


22 


IRAS 19156+1906 


16, 40, 43, 48, 54 


19 17 53.16 


+ 19 12 16.6 


71 


300 


54 


L723 














IRAS 19180+1114 


40, 62 


19 20 25.8 


+ 11 19 52 


71 


300 


3 



L673 A, RNO109 
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Table 3 - continued from previous page 



Object" 


Ref.Obj.'' 


RA (2000) 


Dec (2000) 


Ref.'^ 


d[pc] 


Ref.'^ 


IRAS 19345+0727 


16, 22, 40, 43, 48, 54, 62 


1 q ^7 01 0^ 


-1-07 ^zl 1 q 


71 






B335, LDN663 














IRAS 20050+2720 


40, 49 


90 07 06 1 


1 97 98 ZQ 

-\-A I AO Jy 


ziq 


/ WW 


AQ 
+V 


S 106-SMM 


9, 43 


90 97 1 1 G 
ZU At 11. o 


-1-77 99 A ^ 
1 AA 4j 


1 7 
i / 


AOO 
OWW 


/I 7 


IRAS 20386+6751 


40, 43, 48, 49, 54, 64 


90 "^Q Vif\ ^ 
zU jy UO.J 


1 AO 09 1 'T. 

-\-0o UZ I J 


zi8 

+0 


zLzLO 
'4'4W 




L1157 














GF9-2 


10, 43 


90 ^ 1 "^0 1 
ZU J i J\). i 


-|-uU io jy 


7 1 


900 
ZWW 


A7 


IRAS 2 1017+6742 


48, 62 


91 09 9"^ 1 
Zi WZ Z.J. i 


1 A7 ^/i 1 q 


+o 


Zoo 


A8 
+o 


L1172 














CB230 


68 


91 17 zLO 
All! ^yj.yj 


1 A8 1 7 79 

i^Oo i / J A 


AS 


^JW 


AS 
Oo 


IRAS 2 11 69+6804, LDN 11 77 












L944-SMM1 


62 


91 17 /10 7 
All! +U. / 


-L/l 7 1 8 08 ^ 
+4 J io Uo.J 


A9 

OZ 


700 
/WW 


1 
i 


IC 1396W 


64 


9 1 9A OA 1 
zi ZO UO. i 


-l-<;7 ^A 1 7 
+J / JO i / 


71 
/ i 


7^0 


A/1 
0+ 


IRAS 21246+5743 














CB232 


39 


9 1 ^^7 1 O Q 
zi J / iU.O 


-L/l 7 90 7q A 

++J zu jy.+ 


7q 

iy 


jOW 


97 

A 1 


IRAS 21352+4307 














NGC 7 129 FIRS 2 


33,40 


2143 01.6 


+66 03 26 


50 


1250 


33 


L1251B 


68 


22 38 42.52 


+75 11 45.6 


8 


300 


8 


IRAS 22376+7455 














L1211 


64 


22 47 17.2 


+62 02 34 


42 


725 


64 


IRAS 22451+6154 














CepE-MM 


43, 49, 64 


23 03 13.1 


+61 42 26 


51 


730 


43 


IRAS 23011+6126 














NGC 7538 S 


66 


23 13 44.98 


+61 26 49.2 


66 


3500 


46 


L 1246-SMM 1 


62 


23 25 04.7 


+63 36 40 


62 


730 


62 


IRAS 23238+7401 


48,62 


23 25 46.4 


+74 17 38 


48 


180 


48 


CB244, LDN 12624 














IRAS 23385+6053 


35,43 


23 40 54.5 


+61 10 28 


67 


4900 


43 



" Object name used in this paper. Further common names are given below. 

^ References where the object was mentioned as candidate or confirmed Class source (in the samples listed in Sect. 2). 
In brackets deviations from the general selection criteria are shown. In case the review article from Andre et al. (2000) 
provides the only reference, we additionally list the discovery paper of the source. 

References for the origin of the position (Col. 5) and distance (Col. 7). 
References: (1) Dame & Thaddeus (1985), (2) Clark (1991), (3) Ladd et al. (1991a), (4) Carballo et al. (1992), (5) 
Mezgeretal. (1992a), (6) Prusti et al. (1992a), (7) Zinnecker et al. (1992), (8) Kun & Prusti (1993), (9) Richer et 
al. (1993), (10) Glisten (1994), (11) Persi et al. (1994), (12) Bontemps et al. (1995), (13) Chen et al. (1995a), (14) 
Ward-Thompson et al. (1995b), (15) Wilner et al. (1995), (16) Bontemps et al. (1996), (17) Holland et al. (1996), (18) 
Hurt et al. (1996), (19) Launhardt et al. (1996), (20) Reipurth et al. (1996), (21) Rodriguez & Reipurth (1996), (22) 
Saraceno et al. (1996), (23) Bourke et al. (1997), (24) Chini et al. (1997b), (25) Chini et al. (1997a), (26) Chen et al. 
(1997), (27) Codella & Muders (1997), (28) Launhardt & Henning (1997), (29) Mardones et al. (1997), (30) Olmi et 
al. (1997), (31) Cernicharo et al. (1998), (32) Dent et al. (1998), (33) Eiroa et al. (1998), (34) Hunter et al. (1998), 
(35) Molinari et al. (1998), (36) Andre et al. (1999), (37) Devine et al. (1999), (38) Davis et al. (1999), (39) Huard et 
al. (1999), (40) Park et al. (1999), (41) Reipurth et al. (1999), (42) Tafalla et al. (1999), (43) Andre et al. (2000), (44) 
Chandler & Richer (2000), (45) Gibb & Little (2000), (46) Kalenskii et al. (2000), (47) Looney et al. (2000), (48) 
Shirley et al. (2000), (49) Chini et al. (2001), (50) Fuente et al. (2001), (51) Giannini et al. (2001), (52) Johnstone et 
al. (2001), (53) Lehtinen et al. (2001), (54) Motte & Andre (2001), (55) Persi et al. (2001), (56) Phillips et al. (2001), 
(57) Voung et al. (2001), (58) Ward-Thompson & Buckley (2001), (59) Anglada & Rodriguez (2002), (60) Hartmann 
(2002), (61) Nisini et al. (2002), (62) Visser et al. (2002), (63) Eisloffel et al. (2003), (64) Froebrich et al. (2003), (65) 
Rodriguez et al. (2003), (66) Sandell et al. (2003), (67) Thompson & Macdonald (2003), (68) Young et al. (2003), 
(69) Wolf-Chase et al. (2003), (70) Rengel et al. (2004), (71) SIMBAD database 
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B. Notes to individual objects 



Table 4:: Notes to individual sources 



Object" 



Description 



W3 OH 

L 1448-12 
L1448NW 



L 1448 N 



L 1448 C 
RN015FIR 



RN015 

NGC 133311 
IRAS 03256+3055 



NGC 133312 



SVS 13 B 



NGC 1333 14 A 



The 2MASS source is 3.5" away from the coordinates adopted here. 
Fluxes from O'Linger et al. (1999) taken in the whole observed field are not used. 
L 1448 NW is situated about 20" NNW of the binary L 1448N(A/B). The flux mea- 
surements are strongly influenced by the extended envelope of this nearby source. 
Fluxes from Chandler & Richer (2000) are given in Jy/beam and hence not used. The 
two 2.7 mm fluxes measured in small (< 5") apertures by Looney et al. (2000) are not 
used either. 

This source is actually a protostellar binary consisting of L 1448 IRAS 3A and 
IRAS 3B. Since IRAS 3 A is about one order of magnitude brighter in the FIR as well 
as in the millimeter range, all the unresolved observations are attributed to IRAS 3A. 
Sub-mm fluxes taken in 40" apertures (Shirley et al. (2000)) are not used due to the 
vicinity of L 1448 NW. 

We chose the values measured in boxes larger than 30", in order not to miss emission 
from the outer regions of the extended envelope. 

In the vicinity (about 1' south-east) of RN015FIR the somewhat older source 
RNO 15 can be found. Both objects possess about the same brightness and hence 
flux measurements in large boxes might be influenced by RNO 15. This applies espe- 
cially for the IRAS data. At 60 and 100 fim there are discrepancies between the IRAS 
and ISO data. We used the ISO data since they are measured in a slightly smaller 
aperture. 

The data at short wavelengths and arround the maximum show a large scatter. This 
might partly be due to the nearby source RNO 15 FIR and also due to the fact that most 
of these measurements are taken before 1990 and suffer from large uncertainties. 
Here the fluxes given by Sandell et al. (2001) are excluded. They are unreUable 
because the source is at the edge of their observed field. 

This source seems to consist of three objects (Young et al. (2003)). The IRAS position 
coincides with the weakest of the sub-mm sources. We excluded the fluxes measured 
in the 120" apertures. There is a discrepancy between the 800 and 850 fim point. We 
did not use the 800 /^tm data from Anglada & Rodriguez (2002) since an unusual /3 > 3 
would be obtained. 

NGC 1333 12 is a double source consisting of IRAS 2 A and 2B with a separation of 
about 30". IRAS 2A is much brighter than IRAS 2B in the sub-mm and hence the 
fluxes given for both objects are attributed to IRAS 2A. 

Here it is particular important to have the aperture sizes in which the photometry was 
obtained, since SVS 13 B is a triple protostellar system. It consists of MMS 1, 2, and 
3. While MMS 1 and 2 have about the same brightness, MMS 3 is weaker. There are 
no data for the single sources at wavelengths shorter than 350 fim. 
This is a double source 14 Al and 14 A2. 14 Al is the brighter of the two sources and 
all the photometric data was attributed to this object. 



Continued on next page 
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Table 4 - continued from previous page 



Object" 



Description 



NGC 1333I4B 



IRAS 03282+3035 

HH211MM 
B213 



L1551-IRS5 
L1551-NE 
IRAS 04325+2402 



L1527 

HH114MMS 
IRAS 05173-0555 



LI 634 IRS 7 
RN0 43-MM 

OMC 3-MM6 



OMC 3-MM9 

OMC 3-MM7 

L1641N 
HHl-2 MMS2 

L 1641 VLA 1 



NGC 1333 14 B is situated about 30" south-east of NGC 1333 14 A and hence the 
fluxes measured in large areas might be influenced. The data from Sandell et al. 
(2001) were taken in a relatively small aperture and might not include all the emission 
from the envelope. It was therefore not used. The flux from Rengel et al. (2004) is 
also excluded since the object is at the border of their map. The large fluxes given 
from Rebull et al. (2003) are not used since they might be influenced by surrounding 
emission. 

Even though this source seems to be isolated, the fluxes measured in the 120" box by 
Sandell et al. (2000) are excluded. 

The data from Dent et al. (1998) are taken in an unknown aperture and hence not used. 
We did not use the fluxes measured in a 120" aperture by Young et al. (2003) since 
they might by influenced by other surrounding sources. The IRAS 100 /xm limit was 
used as flux. 

Mesurements taken in very small apertures (< 10") were exluded. Also the data taken 
in the 270" apertures are not used. 

This object is a double source consisting of L 1551-NEA and B. All fluxes are at- 
tributed toL1551-NEA. 

The measurements arround 800 /im and 1.3 mm do not correspond to each other. We 
did not use the low values given by Moriarty-Schieven et al. (1994), since very high 
values of f3 (> 3) are needed to explain these data. Also the upper limits given by 
Barsony & Kenyon (1992) seem to be too low in this context. The 2MASS source is 
6" away from the coordinates adopted here. 

There is a large scatter in the fluxes taken between 350 and 1.3 mm even if they 
are taken with the same aperture size. We do not use data taken in the 120" apertures. 
Difficulties occur since the IRAS points are only upper limits at all wavelengths. The 
1.3 mm flux given by Chini et al. (1997b) for the whole observed field was exluded. 
There is a large scatter in the data at all wavelengths. We did not use the flux of the 
whole field given by Chini et al. (1997b). There are two sets of fluxes in the sub-mm 
data; low and high values. We chose the lower fluxes since they seem to present the 
better the central source and its envelope. 

The 2MASS source is 3" away from the coordinates adopted here. 
The measurement from Chini et al. (1997b) at 1.3 mm for the whole field was ex- 
cluded. 

No data shortward of 100 /im are available. The 185 /um point taken in a 3' aperture 
was not used. All the fluxes below 200 fim seem to be anomolous due to the large 
apertures. 

There is no IRAS detection. The 2MASS source is 6.6" away from the coordinates 
adopted here. 

The upper limit at 100 fim in the IRAS catalogue seems to be too small compared to 

the other data. The 2MASS source is 3.5" away from the coordinates adopted here. 

The 2MASS source is 7" away from the coordinates adopted here. 

There is a large scatter in the 1.3 mm data. We used only the flux measured in the 8.5" 

aperture. 

The scatter in the measurements is very large. We did not use the huge 1.3 mm flux. 
The 2MASS source is 4.3" away from the coordinates adopted here. 
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Table 4 - continued from previous page 



Object" 



Description 



L 1641 SMS III 
HH147MMS 

NGC 2023-MM 1 

NGC2024-FIR5 
HH212-MM 

HH25 MMS 

HH24MMS 

NGC 2068 LBS 17 
HH 1 1 1 MMS 

MonOBlIRAS12Sl 
IRAS 08076-3556 
BHR71-MM 
IRAS 12553-7651 
IRAS 15398-3359 

IRAS 16017-3936 
VLA 1623 



IRAS 16293-2422 

IRAS 16493-4242 
IRAS 18148-0440 

SerpSMM5 
SerpSMM3 
SerpSMM2 

IRAS 18331-0035 

HH108 MMS 
G34.24h-0.13MM 
IRAS 19156H-1906 



IRAS 19180H-1114 

RCrARl 

IRAS 19345H-0727 



The 2MASS source is 6.5" away from the coordinates adopted here. 
Measurements taken in large apertures (> 30") are excluded. The 2MASS source is 
2.5" away from the coordinates adopted here. 

The 450 /um point from Sandell et al. (1999) in the small aperture was excluded, as 

well as the data taken in the 180" apertures. 

The huge 1.3 mm flux in the 30" aperture was excluded. 

The 1.3 mm flux from Chini et al. (1997b), measured in the whole observed field, was 
excluded. 

The two fluxes from Gibb & Davis (1998) measured in large apertures at 60 and 
100 ^m are excluded. The 1 .3 mm point is not used. 

All the IRAS points are just upper limits. The 1.3 mm flux from Launhardetal. (1996) 
measured in a 30" aperture was excluded. 
All data shortward of 100 /im are upper limits. 

The data obtained in large apertures (> 30") were excluded as well as the high flux at 

1.3 mm. The 2MASS source is 4" away from the coordinates adopted here. 

The 2MASS source is 4" away from the coordinates adopted here. 

The 2MASS source is 4.7" away from the coordinates adopted here. 

The 2MASS source is 4" away from the coordinates adopted here. 

The 2MASS source is 7" away from the coordinates adopted here. 

The measurements in the 120" apertures are excluded. The 2MASS source is 3" 

away from the coordinates adopted here. 

The 2MASS source is 9" away from the coordinates adopted here. 
This source is a double object consisting of the two sources VLA 1623 A and B, which 
are comparably bright in the millimeter range. The measurements from Pudritz et al. 
(1996) were excluded. 

IRAS 16293-2422 consits of two sources, IRAS 16293 A and B, seperated by about 
4" (Looney et al. (2000)). We used only the data given for both objects together. 
The 2MASS source is 3.5" away from the coordinates adopted here. 
The data taken in the 19 and 120" apertures are excluded. The 2MASS source is 6" 
away from the coordinates adopted here. 

The 2MASS source is 4" away from the coordinates adopted here. 

The flux at 3.4mm measured in the small aperture was not used. 

The upper limits at long wavelengths are taken within small apertures and hence do 

not contradict the fluxes. 

The 1.3 mm flux taken in a 11" aperture was excluded. The 2MASS source is 3" 

away from the coordinates adopted here. 

The 1.3 mm flux taken in a 1 1" aperture was excluded. 

The 2MASS source is 5" away from the coordinates adopted here. 

The large flux at 1.3 mm given by Cabrit & Andre (1991) was excluded. Also all data 

taken in apertures larger than 100", as well as the other large flux taken at 1.3 mm 

were not used. 

The 2MASS source is 3" away from the coordinates adopted here. 

The 2MASS source is 5" away from the coordinates adopted here. 

IRAS 19345H-0727 shows a large scatter in the sub-mm fluxes. We did not use the data 

taken in apertures larger than 120". 
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Table 4 - continued from previous page 



Object" 



Description 



IRAS 20050+2720 
IRAS 20386+6751 

IRAS 21017+6742 

CB230 
L944-SMM1 
IC1396W 
CB232 



NGC 7 129 FIRS 2 

L1251B 
L1211 



CepE-MM 

NGC 7538 S 
L1246-SMM1 
IRAS 23238+7401 



The 2MASS source is 2.5" away from the coordinates adopted here. 

The sub-mm data taken in apertures larger than 60" were excluded. There is large 

scatter in the sub-mm data. 

Sub-mm fluxes measured in the 120" apertures are excluded. There is a large scatter 
in the sub-mm data. 

The 2MASS source is 2" away from the coordinates adopted here. 

The sub-mm fluxes measured in the 50" apertures by Visser et al. (2002) are used. 

The IRAS data might be confused by other cold dust in the vicinity. 

This is a double source consisting of CB232SMM 1 and SMM2. SMM 1 is much 

brighter in the sub-mm. Hence, all data are attributed to this object. The 2MASS 

source is 6" away from the coordinates adopted here. 

From the 1.3 mm points we chose the 1.49 Jy, since all the others contradict signifi- 
cantly with the data taken at nearby wavelengths. 
Data taken in the 120" apertures are excluded. 

These are several (at least four) sources. Also there are no sub-mm data. The ISO data 
represent the object MMS 4, while the IRAS source might include all objects (MMS 1, 
2, 3, and 4). The ISO data were used and the IRAS points are excluded. The 2MASS 
source is 6" away from the coordinates adopted here. 

This is a double source with a separation of 1.4". The flux of Moro-Martin et al. 

(2001) at 1.3 mm is not used due to the small aperture. 

The 2MASS source is 3" away from the coordinates adopted here. 

The sub-mm fluxes taken in the 50" apertures by Visser et al. (2002) are used. 

The 120" aperture measurments are excluded. 



° Object name used in this paper Further common names are given in Table 3. 
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